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SUMMARY 

A  simple  curve-fitting  method  is  described,  having  application  to 
data  which  obeys  an  assumed  analytic  relationship,  but  in  which  an  unknown 
error  or  offset  is  present  in  the  independent  variable.  The  use  of  the 


method  is  illustrated  by  showing  how  it  may  be  used  in  the  calibration  of 


accelerometers  and  the  analysis  of  drag  measurements. 
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1  INTRODUCTION 

Curve-fitting  procedures  are  widely  used  in  the  analysis  of  experimental  data.  The 

range  of  techniques  employed  is  varied  and  includes  linear  regression  using  least-squares 

1  2  3 

analysis  ,  cubic-spline  fitting  and  various  polynomial-fitting  methods  .  Many  computer 

programs  have  been  written  employing  the  various  curve-fitting  methods  available,  and 

.  4 

some  of  these  have  been  developed  to  a  considerable  degree  of  sophistication  . 

2 

Some  methods,  such  as  cubic-spline  fitting  and  the  method  of  Ref  4,  calculate 

a  curve  which  passes  through  all  the  co-ordinates  supplied;  other  methods  calculate  the 

13 

curve  of  a  specified  form  which  represents  the  'best'  approximation  to  the  data  *  .  The 
procedure  described  in  this  Memorandum  falls  into  the  latter  class  of  curve-fitting 
methods . 

The  method  described  in  this  Memorandum  is  applicable  in  situations  where  a  known 
analytic  relationship  between  the  dependent  and  the  independent  variable  may  be  assumed, 
but  there  exists  some  unknown  offset  error  in  the  measurement  of  the  independent 
variable.  Expressed  in  these  terms  the  method  appears  to  be  rather  limited  in  scope; 
however,  useful  applications  do  exist  and  two  examples  are  discussed  in  this  Memorandum. 

The  method  of  curve  fitting  used  is  discussed  in  section  2  of  this  Memorandum, 
and  examples  of  applications  of  computer  programs  using  this  method  are  discussed  in 
section  3. 

2  CURVE-FITTING  METHOD 

The  method  is  used  to  fit  curves  to  experimental  data  for  which  an  analytic 
relationship  between  the  dependent  and  the  independent  variables  may  be  assumed,  but  in 
which  there  exists  an  unknown  offset  in  the  independent  variable.  That  is  curves  having 
the  form 

y  «  a  f(x  -  e)  +  b  , 


where  e  is  the  unknown  offset  error,  and  a  and  b  are  unknown  constants.  The 
function  f  is  the  analytic  function  assumed  to  apply  to  the  data. 

The  procedure  used  in  the  curve-fitting  method  is  to  use  an  iterative  process  to 
find  e  ,  using  linear  regression  to  control  the  iterations  and  to  determine  the 
coefficients  a  and  b  .  The  iteration  for  c  is  performed  as  part  of  a  computer 
program,  the  method  adopted  being  outlined  below. 

(1)  The  initial  data  is  read.  This  consists  of  the  number  of  co-ordinates  to  be 

input  and  the  range  in  which  e  is  expected  to  lie,  defined  by  the  lowest  and  highest 

values,  e  .  and  e 

mm  max 

(2)  The  co-ordinate  array  of  the  independent  variable,  x  ,  and  the  dependent 
variable,  y  ,  is  stored. 

(3)  An  array  is  created  containing  f(x  -  e  .  )  ,  where  c  .  is  the  initial 

min  mm 

value  of  the  unknown  offset  and  f  is  the  assumed  analytic  function. 


(4)  Linear  regression  analysis  is  performed  on  the  f(x  -  e  ..  )  and  y  arrays, 
yielding  values  of  the  coefficients  a  and  b  ,  and  the  correlation  coefficient  of  the 
straight  line  fitted  to  the  data.  The  equations  used  are  as  follows* 


Z<XY)  -  —g— 

^X2  _  jjy£ 

EY-aEx 


X  and  Y  are  corresponding  values  from  the  f(x  -  e  and  y  arrays  and  N  is 

the  number  of  values  stored  in  the  arrays.  The  correlation  coefficient  is  given  by  the 


equation 


n£(xy)  -£x 

£s>2  -  (rx)2]^2  -(r*)2]1 


(5)  A  new  array  containing  f(x  -  e^)  is  created,  where  is  a  new  value 

for  e  ,  obtained  by  the  equation 


e  =*  t  .  + 

n  min 


(e  -e.)(n-l) 
max  min 


where  n  is  the  iteration  number. 

(6)  Linear  regression  is  again  performed  and  the  correlation  coefficient 
compared  with  that  obtained  in  the  previous  iteration.  If  the  correlation  coefficient 
has  increased  the  iterative  process  is  continued.  Steps  5  and  6  are  repeated  until 

a  correlation  coefficient  is  calculated  which  is  less  than  that  obtained  in  the  previous 
iteration,  or  until  =■  ,  in  which  case  a  warning  is  output  indicating  that  the 

value  of  e  lies  outside  the  specified  range. 

If  the  correlation  coefficient  is  less  than  that  calculated  in  the  previous  itera¬ 
tion  it  is  presumed  that  the  optimum  value  of  e  lies  between  e  and  e  .  .  If  this 

n  n-2 

condition  is  encountered  at  the  second  iteration  the  warning  message  is  output  indicating 
that  e  is  outside  the  specified  range,  and  the  calculation  stops. 

(7)  The  difference  en  -  en_2  is  evaluated.  If  this  is  greater  than  a  prescribed 

amount  e  .  is  set  to  e  .  ,  and  e  is  set  to  e  .  Steps  5  to  7  are  repeated 
min  n-Z  max  n  r 

until  e  -  e  .  is  less  than  the  prescribed  value.  When  this  condition  is  satisfied 
e  is  considered  to  have  been  established  with  sufficient  accuracy.  The  values  of  a, 
b  and  e  are  printed,  together  with  the  correlation  coefficient  and  the  root  mean 
square  error  of  the  curve  fit. 


i  *0  I 
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3  EXAMPLES  OF  APPLICATIONS  OF  THE  METHOD 

The  procedure  set  out  in  the  preceding  section  can  be  applier  with  considerable 
flexibility  since  the  iterative  process  employed  is  independent  of  the  assumed  analytic 
function,  f  . 

Two  examples  of  applications  of  the  curve-fitting  method  are  presented  in  this 
section.  The  first,  discussed  in  section  3.1,  is  concerned  with  the  calibration  of 
accelerometers,  which  are  used  to  determine  the  model  incidence  in  the  5m  wind  tunnel^. 
The  second  example,  described  in  section  3.2,  shows  how  experimental  lift  and  drag  data 
may  be  analysed  to  determine  the  lift-dependent  drag  parameters  of  a  high-lift  wing. 

This  sub-section  also  shows  some  examples  of  the  use  of  the  derived  drag  parameters  in 
data  correction  and  analysis. 

3. 1  Accelerometer  calibration 

The  5m  wind  tunnel*’  is  a  low  speed  pressurised  facility  in  which  large  models 
fitted  with  high-lift  devices  are  tested.  The  loads  generated  are  such  that  significant 
deflections  of  the  models  and  their  support  rigs  may  occur  and  these  deflections  must  be 
accounted  for  in  the  measurement  of  the  model  incidence.  The  method  of  measurement 
adopted  is  to  use  accelerometers  mounted  within  the  model  as  inclinometers  so  that  a 
continuous  measurement  of  the  achieved  model  attitude  may  be  obtained.  In  the  5m  wind 
tunnel  Sundstrand  Q-flex  servo  accelerometers  are  used.  This  type  of  accelerometer  pro¬ 
duces  an  output  voltage  which  is  proportional  to  the  acceleration  component  along  the 
sensitive  axis  of  the  accelerometer. 

Three  accelerometers  are  mounted  in  the  model  in  a  tri-axial  mount,  with  their 
sensitive  axes  nominally  orthogonal  to  each  other^.  The  tri-axial  mount  is  located  in 
a  position  determined  by  the  geometry  of  the  model;  the  accelerometer  axes  and  the  model 
axes  do  not  necessarily  coincide.  For  tests  made  with  the  wings  level,  two  accelerometer 
outputs  are  significant;  those  from  the  accelerometers  which,  when  the  tri-axial  mount  is 
at  zero  incidence,  have  their  sensitive  axes  nominally,  vertical ,  and  nominally  parallel 
with  the  tunnel  axis. 

The  outputs  from  these  two  accelerometers  are  of  the  form 


m.  cos(b  -  s.)  cos(a  -  s.)  +  c, 
i  $  11 


V2  -  m2  cos($  -  s^)  sin(a  -  s2>  +  c2  •  (2 

In  these  equations  V  is  the  output  voltage,  $  the  roll  angle  of  the  model  and  a  is 

the  incidence  of  the  model  datum  plane,  m  is  the  accelerometer  sensitivity,  s,  the 

♦ 

angular  offset  in  roll,  s  the  angular  offset  in  pitch  and  c  is  the  voltage  offset 
(ie  the  output  when  $  ■  s^  and  a  *  s) .  If  the  tri-axial  mount  is  truly  orthogonal 
end  the  accelerometers  have  no  exis  alignment  errors,  Sj  end  s2  will  be  identical. 

In  practice,  however,  the  values  of  s  calculated  in  the  calibration  of  the  two 


accelerometers  are  rarely  identical,  the  error  being  due  partly  to  mounting  and  alignment 
errors  and  partly  to  the  sensitivity  of  the  curve-fitting  method  to  scatter  in  the 
calibration  data. 


In  calibrating  the  accelerometers  it  is  assumed  that 

(a)  tests  are  made  with  the  wings  of  the  model  level.  Thus  $  •  0°  when  the  model 
is  upright  and  180°  when  it  is  inverted; 

(b)  the  angular  offset  in  roll,  s^  ,  is  small. 

With  these  assumptions  the  equations  become 


Vj  =  mj  cos(a  -  s^)  + 

V2  ■=  m2  sin(a  -  s2>  +  c2 
for  tests  with  the  model  upright,  and 


(3) 

(4) 


cos(a  -  Sj)  +  Cj 

(5) 

sin(ct  -  s2)  +  c2 

(6) 

for  tests  with  the  model  inverted. 

The  constants  m,  c  and  s  are  established  for  each  model  test  series  using  computer 

programs  based  on  the  curve-fitting  method  outlined  in  section  2  of  this  Memorandum. 

Calibration  data  for  the  accelerometers  is  obtained  by  recording  their  output  voltages, 

Vj  and  V2  ,  for  a  number  of  incidences  covering  the  anticipated  experimental  range, 

using  an  inclinometer  as  an  incidence  standard.  The  inclinometer  measurements  are 

expected  to  be  accurate  to  1  minute  of  arc,  ie  a  ■  a  ,  ±  0.017°  . 

r  measured 

Fig  1  shows  the  results  of  a  typical  static  calibration  of  the  accelerometers 
fitted  to  a  model  tested  in  the  5m  wind  tunnel.  Fig  la  shows  the  calibration  calculated 
for  the  accelerometer  with  output  equation  (3).  The  calibration  coefficients  obtained 
were  m(  -  3.0516  ,  Sj  -  -2.3184°  and  cf  ■  0.0089  volt.  The  calibration  of  the 
accelerometer  with  output  equation  (4)  is  shown  in  Fig  lb.  The  calculated  coefficients 
for  this  accelerometer  were  m2  -  2.9322  ,  s2  -  -2.3686°  and  c2  -  -0.0076  volt. 

The  fit  of  the  calculated  calibration  curves  to  the  data  was  good,  the  correlation 
coefficients  and  root  mean  square  error  calculated  for  the  two  calibrations  being 


CCj  «  0.9999994  rms  error ^  ■  0.00011  volt 


CC2  “  0.9999988  rms  error2  ■  0.00072  volt. 

For  each  accelerometer  the  significence  of  the  root  mean  square  error  in  volts  varies 
with  incidence.  For  incidences  in  the  range  -45°  <  o  <  45°  the  sine-curve  accelerometer 
(output  equation  (4))  is  sure  sensitive  than  the  cosine-curve  accelerometer  (output 
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equation  (3)),  and  thus  gives  more  accurate  results.  The  root  mean  square  errors  in 
degrees  of  the  two  calibration  curves  are  rms  error ^  -  0.120°,  and  rms  erro^  -0.014°. 
Thus,  despite  a  higher  rms  error  in  volts,  the  sine-curve  accelerometer  achieves  greater 
accuracy  than  the  cosine-curve  accelerometer,  as  expected  since  the  calibrated  range 
was  -6°  <  a  <  30°  . 

When  a  model  is  to  be  tested  in  both  the  upright  and  the  inverted  configurations 
the  same  values  of  m,  c  and  s  should  be  used  for  both  cases,  since  the  calibration 
coefficients  are  unaffected  by  inversion  of  the  model.  However,  independent  evaluation 
of  the  calibration  coefficients  will  not  necessarily  result  in  an  identical  set  being 
calculated  for  the  upright  and  inverted  cases,  due  to  the  inevitable  scatter  present  in 
the  calibration  data. 

The  introduction  of  the  roll  angle,  ,  into  the  accelerometer  output  equations 

t  ... 

allows  a  single  set  of  accelerometer  output  constants  to  be  derived  from  data  obtained  in 
separate  calibrations  witn  the  model  in  the  fipright  and  the  inverted  configurations. 
Output  equations  (1)  and  (2)  are  used,  with  assumed  negligible  and  $  set  to  0°  for 

the  upright  data  and  180°  for  the  inverted  data. 

Fig  2  shows  the  results  of  such  a  calibration  for  a  sine-curve  accelerometer 
( ie  obeying  output  equation  (2)}.  The  calibration  data  was  obtained  for  an  irtcidence 
range  of  -6°  <  a  <  6°  for  both  the  upright  and  inverted  cases.  The  output  equation 
calculated  from  the  calibration  data  was 


3 


V  «  2.9692  cos  <t>  sin(a  -  0.0086°)  +  0.00666  , 


and  the  rms  error  of  the  calibration  curve  was  0.0071°.  The 
in  gradient  of  the  output  curve  associated  with  inversion  of 
a  positive  gradient  was  obtained  with  the  model  upright  (*  = 
gradient  was  obtained  with  the  model  inverted  ($  =  180°). 


figure  illustrates  the  change 
the  model;  the  data  having 
0°) ,  and  that  with  a  negative 


3.2  Lift-dependent  drag  analysis 

The  drag  polar  of  an  uncambered,  untwisted  wing  may  be  represented  by  the  simple 
drag  equation 

c  =  c  +  — 

CD  CDq  nA  °L  * 


where  is  the  profile  drag  of  the  wing,  which  is  the  minimum  drag  and  occurs  at 

zero  lift  if  the  simple  drag  equation  applies.  The  term  (K/uA)C^  is  the  drag  due  to 
lift,  K  being  the  lift-dependent  drag  factor,  and  A  the  aspect  ratio  of  the  wing. 

This  drag  equation  is  an  idealised  representation  of  the  drag  polar  and  is  applic¬ 
able  only  to  a  limited  range  of  lift  coefficients.  At  high  lift  coefficients,  as 
Cf  ^  is  approached,  the  occurrence  of  flow  separations  from  the  surface  of  the  wing 
causes  the  drag  to  increase,  invalidating  this  simple  equation.  When  the  lift  coefficient 
is  snail  the  laminar  region  of  the  boundary  layer  on  the  wing's  surface  may  have  a 
significant  extent,  causing  a  reduction  in  the  drag  and  again  invalidating  the  simple 
drag  equation. 


*  j 
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If  the  wing  is  cambered  the  drag  equation  must  be  modified  since,  although  the  drag 

polar  remains  parabolic,  the  minimum  drag  no  longer  occurs  at  zero  lift,  but  at  some  other 

lift  coefficient,  Cr  .  .  The  modified  form  of  the  equation  is 

’  '-min 


The  applicability  of  this  equation  is  affected  by  the  same  limitations  as  the  simple 
drag  equation  and,  additionally,  if  the  cambered  wing  is  fitted  with  high-lift  devices, 
further  limitations  may  be  imposed  by  the  occurrence  of  flow  separations  from  these 
devices . 

A  computer  program  has  been  written  to  calculate  the  drag  parameters  Cn  .  »  K 

and  Ct  .  from  measured  lift  and  drag  coefficients,  using  the  calculation  method  set 
tmin 

out  in  section  2  of  this  Memorandum.  The  program  depends  on  the  hypothesis  that  the 
modified  drag  equation  set  out  above  is  valid  for  at  least  a  portion  of  the  measured  drag 
polar,  and  that  sufficient  data  points  have  been  obtained  in  the  region  of  validity  to 
allow  a  satisfactory  fit  to  the  data  to  be  obtained. 

The  drag  parameters  are  important  because,  once  calculated,  they  determine  the 
attached-flow  drag  polar,  which  serves  as  a  general  description  of  the  aerodynamic 
performance  of  the  aircraft  or  model  for  which  it  is  established,  in  the  range  of 
validity  of  the  drag  equation  used.  The  attached-flow  drag  polar  also  has  certain 
specific  applications,  as  noted  below. 

(a)  The  initial  occurrence  of  flow  separation  and  the  changes  in  drag  consequent  on  the 
flow  separation  may  be  established  by  comparison  of  the  attached-flow  and  the  measured 
drag  polars.  The  changes  in  drag  observed  may  be  increments,  if  the  separation  occurs 

at  high  incidence,  or  decrements  if  the  separation  is  from  the  rear  surface  of  a  slat  at 
low  incidence. 

(b)  In  wind  tunnel  testing,  knowledge  of  the  separated  drag  is  required  if  the  blockage 
effects  due  to  the  separated  wake  are  to  be  accounted  for  using  the  method  recommended 

by  Maskell  .  Experience  has  shown  that  to  make  an  objective  assessment  of  the  separated 
drag  for  a  model  with  high-lift  devices  deployed  the  attached-flow  drag  polar  must  be 
represented  by  the  modified  drag  equation  quoted  previously.  The  calculation  method 
described  in  this  Memorandum  was  initially  developed  for  this  application. 

(c)  Comparison  of  attached-flow  drag  polars  may  be  useful  in  the  examination  of  the 
effects  of  changes  in  high-lift  device  deployment  angles,  and  the  assessment  of  incremental 
effects  such  as  those  due  to  the  carriage  of  stores,  the  lowering  of  the  undercarriage  of 
an  aircraft,  or  the  fitting  of  gauzes  to  restrict  the  flow  through  the  engine  nacelles  of 

a  wind  tunnel  model. 

Fig  3  shows  the  lift  and  drag  curves  of  a  model  of  a  commercial  transport  aircraft 
tested  in  the  5m  wind  tunnel  in  a  take-off  configuration,  with  full-span  leading  edge  slats 
deployed  and  the  flap  system  retracted.  Drag  parameters  were  derived  from  the  data  shown 
in  Fig  3  and  used  to  define  the  attached-flow  drag  polar.  The  drag  parameters  were 
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determined  by  curve-fitting  the  lift  and  drag  data  in  the  range  °-‘  <  <VW '  °-8  • 
The  attached-flow  drag  polar  is  defined  by  the  equation 


CD  '  CD  .  *  S  °L  -  CL  . 

un  \  min/ 

where,  for  this  set  of  data,  0^^  =  0.0427  ,  =  0.2291  ,  K  =  1.4951  and 

A  ■  7.73  .  The  measured  polar  and  the  calculated  attached-flow  drag  polar  are  shown 
in  Fig  4. 

Comparison  of  the  attached-flow  polar  and  the  measured  data  shows  that  the  initial 
onset  of  flow  separation  occurred  at  about  CL  =  1.5  .  At  lift  coefficients  of  1.62  and 
1.64  increments  in  drag  reveal  the  occurrence  of  flow  separations  from  the  two  underwing 
engine  nacelles  fitted  to  the  model,  and  the  separated  drag  then  increases  rapidly  as  the 
lift  coefficient  approaches  • 

At  lift  coefficients  between  C  =  -0.1  and  C  =  0.6  the  drag  of  the  model  is 

Ij  la 

less  than  that  suggested  by  the  attached-flow  drag  polar.  At  the  low  angles  of  incidence 
associated  with  this  lift  coefficient  range'  (-4°  <  a  <  +1°),  flow  separation  occurs  from 
the  rear  surface  of  the  slat,  reducing  the  pressure  drag  which  would  otherwise  have 
arisen  due  to  the  suctions  on  this  surface.  A  further  effect  of  this  flow  separation  is 
a  slight  increase  in  the  lift  coefficient  at  low  angles  of  incidence  compared  to  that 
anticipated  from  the  linear  portion  of  the  lift  curve.  This  effect,  which  may  be  seen  in 
Fig  3,  is  due  to  the  reduced  downward  component  of  the  slat  load  in  the  presence  of  the 
flow  separation. 

Ultimately,  for  angles  of  incidence  less  than  -4°,  the  flow  separation  from  the 
rear  surface  of  the  slat  develops  to  such  an  extent  that  the  increase  in  drag  due  to  the 
separated  wake  outweighs  the  reduction  in  the  pressure  drag,  resulting  in  the  net  increase 
in  drag  observed  at  negative  lift  coefficient  in  Fig  4. 

In  order  to  correct  wind  tunnel  measurements  for  the  effects  of  the  blockage  due 
to  the  wake  of  the  model  and  its  support  rig,  it  is  necessary  to  be  able  to  derive  the 
lift-dependent  drag  parameters  from  the  measured  data.  The  correction  due  to  the  wake 
blockage  may  be  calculated  using  an  equation  of  the  form^ 


(\  *  \)  *  “  (\J  ' 


In  this  equation  6(q/qQ)  is  the  correction  to  be  applied  to  the  wind  tunnel  dynamic 

pressure,  S  is  the  model  wing  area  and  C  the  tunnel  cross-sectional  area.  CpQ  is 

the  profile  drag  of  the  model  and  occurs  at  -  CLq  .  The  profile  drag  is  assumed  to 

be  the  minimum  drag  in  attached  flow,  and  hence  Cn,,  *  Co  .  and  Cr „  -  Ct  .  .  Cd_  i 

u  "min  u  min  1 

the  tare  and  interference  drag  of  the  model  supports  and  Cd#  i-8  the  separated  drag 
coefficient,  given  by 

CD  CD  "  CD„  “  °D  "  CDI 

sep  T  0  I 


I 


Q> 
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where  CDj.  is  the  lift-dependent  drag 


Fig  5  shows  the  variation  of  the  measured  drag  with  (Cl  -  CLq)  for  the  data  of 
Fig  2,  and  illustrates  the  definitions  of  the  parameters  CDq  ,  and  Cds  '  Als0 

shown  in  Fig  5  is  the  calculated  attached-flow  drag  polar,  which  is  seen  to  be  a  very  good 

fit  to  the  data  in  the  range  for  which  the  curve-fitting  process  was  carried  out, 

2 

0.2  <  (Cl  -  Clq)  <  1.3  .  The  root  mean  square  error  of  the  attached-flow  drag  polar 

in  this  range  is  AC_,  .  =  0.0002  . 

D(rms) 

A  further  application  of  the  curve-fitting  program  to  the  analysis  of  drag  data  is 
illustrated  in  Fig  6,  which  shows  the  effect  of  stores  carriage  on  the  drag  of  a  mode1  of 

g 

a  ground  attack/ training  aircraft  .  Data  is  presented  in  the  Figure  for  four 
conf igurat ions : 

(a)  the  clean  configuration, 

(b)  with  two  stores  pylons  carried  on  each  wing, 

(c)  with  droptanks  on  the  inboard  pylons  and  the  outboard  pylons  fitted  but 
carrying  no  stores,  and 

(d)  with  droptanks  on  the  inboard  pylons,  and  two  rocket  packs  fitted  to  a  twin 
store  carrier  on  each  outboard  pylon. 

For  each  of  these  configurations  the  measured  drag  coefficient  is  plotted  against 
2 

(C.  -  Cl  .  )  using  the  value  of  Cl  •  calculated  for  each  set  of  data  using  the 
L  ‘'min  °  ‘'min  ° 

curve-fitting  program.  The  lift-dependent  drag  parameters,  defining  the  attached-flow 
drag  polar  for  each  run,  are  also  shown  in  the  Figure. 

The  increments  in  the  profile  drag  coefficient  and  the  lift-dependent  drag  factor 
resulting  from  the  carriage  of  the  stores  are  established  in  an  objective  fashion  by  the 
curve-fitting  program,  making  the  technique  particularly  useful  for  this  type  of  com¬ 
parative  analysis. 

4  CONCLUSIONS 

The  curve-fitting  method  described  in  this  Memorandum  forms  the  basis  of  computer 
programs  used  to  calibrate  accelerometers  and  to  derive  attached-flow  drag  polars  from 
measured  lift  and  drag  coefficients.  The  results  presented  in  this  Memorandum  show  that 
the  curve-fitting  method  works  well  and  illustrate  some  of  the  uses  of  the  data  derived 
using  the  drag  analysis  computer  program.  Details  of  the  various  computer  programs 
based  on  this  curve-fitting  method  may  be  obtained  through  Aerodynamics  Department, 

RAE,  Farnborough. 
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Fig  4  Measured  and  attached-flow  drag  polars 
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Fig  5  Parameters  used  in  wake  blockage  calculation 
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